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We  adapt  a  general-purpose  optimization-based  parameter  estimation  technique  previously  described 
in  the  literature  [1  ]  to  evaluate  the  suitability  of  a  number  of  common  kinetic  models  for  the  representa¬ 
tion  of  key  performance  characteristics  (conversion  and  selectivity)  of  catalysts  used  for  the  preferential 
oxidation  of  CO  in  the  presence  of  H2.  We  find  that,  for  process  engineering  applications,  there  is  no  clear 
practical  advantage  to  using  mechanistically  based  kinetic  models  (e.g.  Langmuir-Hinshelwood)  unless 
the  precise  chemical  mechanism  is  known.  Empirical  rate  models  are  found  generally  to  provide  equiv¬ 
alent  or  better  simulations  of  key  performance  variables  for  a  diverse  group  of  catalyst  formulations. 
Furthermore,  we  demonstrate  that  the  water-gas-shift  (WGS)  reaction  is  relevant  within  PROX  reaction 
systems  under  conditons  containing  high  fractions  of  C02  and  H2,  confirming  the  expectations  of  Choi 
and  Stenger  (2004)  [2].  Finally,  we  attempt  to  identify  any  emergent  trends  in  kinetic  parameters  among 
catalysts  sharing  similar  active  metal  or  metal  oxide  components.  Unfortunately,  apart  from  confirming 
that  the  activation  barrier  for  CO  oxidation  is  generally  less  than  the  barrier  for  H2  oxidation  (an  expected 
relationship  for  PROX  catalysts),  no  such  trends  are  found. 

©  2012  Published  by  Elsevier  B.V. 


1.  Introduction 

Noble  metals  (platinum  and  palladium,  in  particular)  and  noble 
metal  alloys  (platinum/ruthenium  blends)  remain  the  anode  cat¬ 
alyst  of  choice  in  state-of-the-art  polymer-electrolyte-membrane 
(PEM)  fuel-cell  devices.  This  is  because  no  other  material  has  yet 
been  shown  to  provide  equivalent  electrochemical  performance 
[3].  Unfortunately,  noble  metal  catalysts  are  well  known  to  be 
susceptible  to  poisoning  by  both  sulfur  and  carbon  monoxide,  espe¬ 
cially  at  lower  temperatures  (<100°C).  Therefore,  to  both  extend 
the  useful  life  of  the  fuel-cell  device  and  to  reduce  the  materials 
cost  associated  with  replacing  expensive  noble  metal  electrodes, 
it  is  crucial  ensure  that  catalyst  poisoning  is  minimized.  The  exact 
amount  of  carbon  monoxide  that  can  be  tolerated  for  fuel  cell  appli¬ 
cations  varies  with  the  amount  of  precious  metal  used  and  with 
the  operating  temperature  of  the  device,  but  is  typically  less  than 
10  ppm  for  PEM  fuel  cells. 
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It  is  estimated  that  80%  of  the  hydrogen  produced  world¬ 
wide  comes  from  the  reforming  of  methane  and  other  light 
hydrocarbons.  The  hydrogen-rich  gaseous  mixture  produced  by 
such  reforming  processes  typically  contains  between  80,000  and 
120,000  ppm  (8-12%)  carbon  monoxide.  The  carbon  monoxide 
content  is  often  reduced  (and  hydrogen  content  simultaneously 
enriched)  using  a  follow-on  water-gas-shift  reaction,  but  the 
amount  of  carbon  monoxide  in  the  final  product  stream  is  typically 
about  10,000  ppm  (1%). 

The  preferential  (or  selective)  catalytic  oxidation  of  carbon 
monoxide  to  carbon  dioxide  (PROX)  is  considered  to  be  an  eco¬ 
nomical  and  highly  distributable  process  for  the  removal  of  trace 
amounts  of  carbon  monoxide.  A  large  number  of  catalyst  formula¬ 
tions  for  this  process  have  been  described  in  the  literature,  most 
of  which  involve  noble  metals  (platinum,  palladium,  ruthenium,  or 
rhodium)  as  the  active  element.  In  the  case  of  PROX  processes,  noble 
metal  catalyst  poisoning  by  carbon  monoxide  is  less  of  a  concern, 
as  the  reactors  typically  operate  at  moderately  high  temperatures 
(between  150°C  and  250  °C). 

In  the  present  work,  we  leverage  a  robust  general-purpose 
optimization-based  parameter  estimation  technique  previously 
described  in  the  literature  [1]  to  evaluate  the  suitability  of  dif¬ 
ferent  kinetic  models  for  the  accurate  simulation  of  PROX  catalyst 
performance  characteristics:  specifically  the  overall  CO  conversion 


AC.  Baughman  et  al.  /  Journal  of  Power  Sources  210  (2012)  402-408 


403 


Nomenclature 

JZC  the  rate  of  carbon  monoxide  oxidation,  as  modeled 
using  an  empirical  rate  law 

JZh  the  rate  of  hydrogen  oxidation,  as  modeled  using  an 
empirical  rate  law 

JZl  the  rate  of  carbon  monoxide  oxidation,  as  modeled 
using  a  mechanistic  (Langmuir-Hinshelwood-type) 
rate  law 

JZD  the  rate  at  which  the  water-gas-shift  reaction 
occurs 

kc  the  pre-exponential  rate  constant  (parameter) 

appearing  in  JZC 

kH  the  pre-exponential  rate  constant  (parameter) 

appearing  in  JZH 

/<Lii  the  pre-exponential  rate  constant  (parameter) 

appearing  in  the  numerator  of  JZl 
kLi 2  the  pre-exponential  rate  constant  (parameter) 

appearing  in  the  denominator  of  JZi 
kD  the  pre-exponential  rate  constant  (parameter) 

appearing  in  JZd 

JCeq  the  equilibrium  constant  for  the  water-gas-shift 
reaction 

Ec  the  activation  energy  (parameter)  appearing  in  JZC 
Eh  the  activation  energy  (parameter)  appearing  in  JZh 

/q,  i  the  activation  energy  (parameter)  appearing  in  the 

numerator  of  JZl 

kLy 2  the  activation  energy  (parameter)  appearing  in  the 

denominator  of  JZL 

ko  the  activation  energy  (parameter)  appearing  in  JZd 
a  the  order  (parameter)  of  JZC  with  respect  to  carbon 

monoxide 

b  the  order  (parameter)  of  JZC  with  respect  to  oxygen 

c  the  order  (parameter)  of  JZh  with  respect  to  oxygen 

Cco  the  local  concentration  of  carbon  monoxide 
Co2  the  local  concentration  of  oxygen 
Ch2o  the  local  concentration  of  water  vapor 
R  the  universal  gas  constant. 

T  temperature 

A  the  cross-sectional  area  of  the  continuous-flow 

reactor 

V  the  superficial  flow  velocity  of  the  reaction  mixture 

Q  the  space-veloctiy  of  the  reaction  mixture.  Calcu¬ 

lated  as  V/L 

T  the  estimator  function 

V  the  set  of  parameter  values  to  be  estimated 

M  the  number  of  operating  temperatures  for  which 
conversion  or  selectivity  data  is  available 
Xk  the  conversion  achieved  by  the  catalyst  operating  at 

temperature  Tk 

Sk  the  oxidative  selectivity  achieved  by  the  catalyst 
operating  at  temperature  Tk 

Rx  the  coefficient  of  determination  for  the  model  fitting 

of  conversion  data 

Rs  the  coefficient  of  determination  for  the  model  fitting 

of  selectivity  data 


achieved,  and  their  selectivity  towards  CO  oxidation  versus  H2  oxi¬ 
dation.  To  do  so,  we  first  formulate  a  general  dynamic  material 
balance  model  describing  a  gas-phase  isothermal  plug-flow- 
reactor  (PFR)  typical  of  those  commonly  used  in  such  processes. 
Then,  four  different  kinetic  (rate)  models  are  developed,  based  upon 
use  of  an  emprical  or  mechanistic  approach,  with  and  without  the 
simultaneous  participation  of  a  water-gas-shift  reaction.  Finally, 


using  our  parameter  estimation  technique,  each  kinetic  model  is 
fit  to  the  experimental  conversion  and  selectivity  measurements 
taken  from  several  reports  in  the  recent  literature. 

We  close  by  evaluating  the  effectiveness  by  which  each  model 
recaptures  catalyst  performance,  discussing  the  significance  of  the 
water-gas-shift  reaction  within  PROX  processes,  and  attempting 
to  identify  any  evident  kinetic  parameter  trends  according  to  the 
active  catalyst  component. 

2.  Computational  methods 

2.1.  Reactor  model 

Classically,  the  preferential  oxidation  of  CO  considers  the  simul¬ 


taneous  catalysis  of  two  reactions: 

CO  +  io2  ->  C02  AH°298  =  -283.0  kj  mol"1  ,  (1) 

H2  +  lo2  -*  H20  AH°98  = -241.8  kj  mol"1.  (2) 

The  kinetics  of  these  two  reactions  have  been  modeled  in  the 
literature  [2,4,5]  using  generalized  empirical  rate  laws  of  the  form: 

i?c  =  fccexp(-|0q!oc£2,  (3) 

rh  =  kH  exp  (- 1^)  Cq2  .  (4) 


It  has  been  shown  [6-8]  that  the  apparent  rate  of  the  hydrogen 
oxidation  reaction  is  independent  of  the  local  hydrogen  concentra¬ 
tion,  and  hydrogen  is  therefore  omitted  from  Eq.  (4). 

The  catalytic  oxidation  of  carbon  monoxide  is  generally  believed 
to  proceed  through  a  surface-reaction  mechanism,  which  has 
prompted  a  number  of  researchers  to  model  the  kinetics  of  reac¬ 
tion  1  using  more  mechanistically  oriented  rate  laws  [9-11].  One 
such  law,  examined  by  Venderbosch  and  co-workers  [12],  follows 
a  Langmuir-Hinshelwood  (LH)  type  mechanism: 

R  _  exp(-(EL1/RT))CCoCo2  ^ 

i_(l+fci,2  exp(-(EL,2/RT))CCo)2 ' 

Additionally,  it  has  been  suggested  [2]  that  the  water-gas-shift 
(WGS)  reaction  can  occur  under  the  conditions  present  within  a 
typical  PROX  reactor. 

H20  +  CO  ^  C02  +  H2  AH298  =  -41.1  kj  mol"1,  (6) 

More  specifically,  large  amounts  of  hydrogen  and  carbon  dioxide 
commonly  present  in  the  reactant  gas  stream  could  promote  the 
reverse  of  reaction  (6),  wherein  carbon  monoxide  is  regenerated. 
It  has  been  shown  [13]  that  the  kinetics  of  this  process  are  well 
described  using  the  following  rate  law: 

RD  =  kD  exp  (-§)  (cC0Ch20  -  •  (7) 

The  equilibrium  constant  of  the  WGS  reaction  is  known  to  be 
well  estimated  as:  Keq  =  exp  ((4577.8 /T)  -  4.33)  [14]. 

The  typical  configuration  for  a  PROX  reactor  is  a  fixed-bed  tubu¬ 
lar  continuous-flow  system.  For  any  chemical  species  i  reacting 
within  such  a  system,  the  material  balance  for  i  surrounding  an 
3z-thick  axial  segment  of  the  reactor  can  be  constructed  as  follows: 

^(QAdz)  =  RtAdz  +  A(  VCt \z  -  VQ |z+jfe).  (8) 

Q  is  the  local  molar  concentration  of  species  i,  A  is  the  cross- 
sectional  area  of  the  reactor,  Rj  is  the  kinetic  rate  law  describing  the 
consumption  or  generation  of  species  i,  and  V  is  the  local  superficial 
flow  velocity  of  the  reacting  mixture. 
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Table  1 

Expressing  the  rate  at  which  each  individual  species  is  consumed  or  generated,  for  each  kinetic  model,  as  a  stoichiometric  sum  of  reaction  rates.  7 Zc  is  described  in  Eq.  (3), 
Rh  in  Eq.  (4),  RL  in  Eq.  (5),  and  RD  in  Eq.  (7). 


Rate 

Kinetic  model 

Empirical 

Empirical  +  WGS 

Mechanistic 

Mechanistic  +  WGS 

Rco 

-Rc 

—Rc  —  Rd 

-Rl 

—Ri  —  Rd 

Rco2 

Rc 

Rc  +  Rd 

Rl 

Rl  +  Rd 

Rh2 

-Rh 

—Rh  +  Rd 

-Rh 

—Rh  +  Rd 

Rh2o 

Rh 

Rh  -  Rd 

Rh 

Rh  -  Rd 

Ro2 

l/2(— 7^c-7^H) 

\/2(-Rc-Rh) 

\/2(-Rc-Rh) 

1/2  (~rc-Rh) 

Kinetic  parameters 

7 

9 

1 

9 

The  particular  ft,  for  each  reacting  species  within  the  PROX  sys¬ 
tem  depends  upon  the  particular  rate  laws  chosen  to  model  the 
reaction  system.  As  we  have  discussed,  there  is  support  in  the  lit¬ 
erature  for  the  modeling  of  carbon  monoxide  oxidation  using  both 
the  generic  empirical  rate  law  given  in  Eq.  (3),  and  more  mecha¬ 
nistic  rate  laws  such  as  that  given  in  Eq.  (5).  It  is  also  not  clear  to 
what  extent  the  WGS  (or  reverse  WGS)  reaction  is  in  fact  active 
within  the  PROX  reactor.  For  these  reasons,  in  the  present  investi¬ 
gation,  we  chose  to  consider  four  alternate  possible  cases  for  each 
ft,.  These  correspond  to  the  use  of  empirical  or  mechanistic  rate 
laws  describing  carbon  monoxide  oxidation,  each  with  or  without 
simultaneous  participation  of  the  WGS  reaction.  The  appropriate 
ft, ’ s,  expressed  as  the  proper  stoichiometric  sum  of  rate  laws,  are 
shown  in  Table  1. 

Under  steady-state  conditions,  Eq.  (8)  neatly  simplifies  to: 

<9> 

Rigorously,  the  magnitude  of  the  local  flow  velocity  depends 
upon  the  local  temperature,  pressure,  and  molar  density  of  the 
reaction  mixture.  Of  these  three,  temperature  is  generally  of  least 
concern  for  modeling  purposes  since  PROX  reactors  are  customar¬ 
ily  operated  under  isothermal  conditions.  Furthermore,  pressure 
drops  due  to  flow  through  the  catalytic  bed  are  typically  quite 
small  (<5kPa  [5],  <0.1  atm  [11]).  It  is  therefore  not  unreasonable 
to  assume  isobaric  operation.  Finally,  even  though  both  oxidation 
reactions  cause  a  net  decrease  in  the  total  mole  count,  the  overall 
effect  of  these  reactions  upon  molar  density  within  the  reactor  is 
minimal.  This  is  because  both  the  total  amount  of  CO  entering  the 
reactor  («1%),  and  the  overall  extent  of  H2  oxidation  are  usually 
each  quite  small. 

As  a  result,  we  assume  that  the  flow  velocity  within  the  reactor 
may  be  considered  constant. 

Knowing  this,  and  making  a  substitution  in  variables  (q  =  z/L, 
where  L  is  the  length  of  the  tubular  reactor)  we  may  re-express 
Eq.  (9)  as: 


=  R,.  =*.  Al  =  IR 

dq  dq  Q 


(10) 


This  alternative  representation  of  the  material  balance  is  advan¬ 
tageous  for  modeling,  as  it  is  independent  of  specific  reactor 
geometry.  The  variable  q  is  a  dimensionless  fraction  of  overall  reac¬ 
tor  length  (0  <  q  <  1 ),  and  the  quantity  V/L  =  Q.  is  the  space  velocity 
of  gases  within  the  reactor  (i.e.  the  number  of  reactor  volumes 
processed  per  unit  time). 

This  material  balance  makes  additional  implicit  assumptions: 
(1)  that  the  flow  velocity  is  sufficently  large  such  that  convec¬ 
tive  mass  transfer  dominates  any  diffusive  mass  transfer  processes 
along  the  long  axis  of  the  reactor,  (2)  that  this  velocity  is  uniform 
over  the  cross-section  of  the  reactor,  and  (3)  that  flow  is  sufficiently 
turbulent  so  as  to  cause  perfect  mixing  in  the  radial  (non-axial) 
dimension.  In  general,  all  three  of  these  assumptions  are  valid  for 
reactors  operating  under  sufficiently  high  space  velocities  (>1  vol¬ 
ume  per  second). 


2.2.  Basis  for  parameter  estimation 

The  usefulness  of  any  mathematical  model  hinges  upon  its  abil¬ 
ity  to  accurately  represent  real  systems.  For  the  model  we  have  just 
described  to  accurately  represent  real  PROX  systems,  we  are  faced 
with  the  need  to  identify  suitable  values  for  those  parameters  (V) 
appearing  in  the  rate  laws. 

We  previously  described  the  use  of  a  straight-foward  general- 
purpose  optimization-based  technique  for  the  estimation  of 
unknown  parameter  values,  applied  to  dynamic  models  of  cell 
metabolism  [1].  The  technique  applies  numerical  discretization 
techniques  to  convert  differential  equations  comprising  the  model 
of  interest  into  an  approximating  set  of  purely  algebraic  expres¬ 
sions.  The  discretized  expressions  are  then  used  to  constrain  an 
optimization  problem  whose  objective  is  to  minimize  the  error 
between  system  performance  as  predicted  by  the  model  against 
that  observed  experimentally. 

Numerically  discretizing  a  generic  function  #  assumes  that 
when  given  a  particular  function  value  corresponding  to  a  partic¬ 
ular  independent  variable  value  (q),  then  the  function  value  at  a 
different  arbitrary  value  of  the  independent  variable  ( q  +  8q)  can  be 
estimated  using  an  algebraic  relationship  of  the  form: 

*1  ,+s,"  ^,=^(^9*^)-  (11) 

The  nature  of  the  estimator  T  defines  the  chosen  discretiza¬ 
tion  method.  Arguably  the  simplest  discretization  technique  is  the 
explicit  first-order  Euler  method,  wherein  T  =  dft/dq  |  ^ 8q .  Although 
quite  straight-forward,  estimates  generated  using  this  method  can 
occasionally  have  substantial  error. 

From  the  fundamental  theorem  of  the  calculus,  it  is  known  that: 


Comparing  Eq.  (12)  against  Eq.  (11)  it  is  clear  that  a  per¬ 
fect  estimator  T  precisely  calculates  the  area  under  the  curve 
d#/dq  between  q  and  q  +  8q.  When  using  an  explicit  first-order 
Euler  estimator,  the  left-hand  panel  of  Fig.  1  shows  the  poten¬ 
tially  large  difference  (error)  between  the  estimated  and  actual 
values  of  this  area  for  a  generic  d#/dq.  The  right-hand  panel 
of  this  figure,  however,  shows  how  the  calculation  can  be 
improved  using  a  simple  second-order  trapezoidal  estimator: 
T=  l/2(d#/dq|^  +  d#/dq|  q+5(?)&7-  This  approach,  a  second-order 
Runge-Kutta  method,  is  sometimes  called  Heun’s  method. 

Because  of  the  potential  for  reduction  in  estimation  error  offered 
with  only  a  minor  increase  in  computational  complexity,  this 
second-order  method  was  chosen  for  the  present  parameter  esti¬ 
mation. 

2.3.  Parameter  estimation  problem  formulation 

The  performance  of  a  PROX  catalyst  is  typically  measured  by 
the  carbon  monoxide  conversion  (Xk  =  Ccolo,rk  -  Ccoli,Tk/Ccolo,rk) 
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Fig.  1.  Illustrating  the  potential  reduction  in  error  by  using  a  simple  second-order  estimator. 


and  oxidative  selectivity  (Sk  =  l/2((Ccolo,rk  -  Cco  li  ,rfc )/( Co2 10  Tfc  - 
Co-,  L  ^  )))  achieved  at  a  particular  operating  temperature  Tk.  We 

z  1  Mfe 

therefore  wish  our  model  to  reproduce  real  (observed)  conversion 
and  selectivity  values  as  closely  as  possible.  Mathematically,  we 
phrase  this  goal  as  maximizing  the  quantity: 


r2x+r% 


where 

r2  =  1  Si^-4)2 


(13) 


(14) 


and  similarly  for  R2. 

Each  R 2  quantity  is  a  coefficient  of  determination,  measuring 
the  quality  by  which  conversion  and  selectivity  value  estimated 
using  the  kinetic  parameters  in  our  model  ( Xk  and  Sk )  capture  the 
variation  in  the  experimentally  observed  data  {S°k,  S°k),  as  compared 
to  simply  using  the  experimental  mean  (X°k,  S°k)  as  the  estimate.  M 
is  the  number  of  distinct  temperatures  for  which  conversion  and 
selectivity  data  have  been  observed. 

For  any  particular  temperature  Tk,  the  discretized  material  bal¬ 
ance  for  species  i  is  represented  using  the  second-order  estimator 
as  follows: 


parameters  for  PROX  catalysts  described  in  eight  different  investi¬ 
gations  taken  from  the  recent  literature.  The  chosen  investigations 
are  summarized  in  Table  2,  and  cover  a  variety  of  catalyst  for¬ 
mulations  commonly  found  in  practice  (e.g.  supported  platinum, 
copper/copper-oxide,  gold,  and  ruthenium)  [16]. 

All  of  the  investigations  chosen  report  experimentally  observed 
values  for  carbon  monoxide  conversion  as  a  function  of  tempera¬ 
ture  (Ap.  The  number  of  distinct  conversion  values  reported  (and 
thus  the  number  of  distinct  operating  temperatures  considered)  in 
each  investigation  is  listed  in  the  table  {M).  Additionally,  the  inlet 
gas  feed  composition  (i.e.  C?)  and  operating  space  velocity  (Q)  are 
also  listed  for  each  investigation. 

Not  every  investigation  we  considered  provided  values  for  the 
oxidative  selectivity  of  their  catalyst  (Sk).  These  investigations 
are  indicated  in  Table  2.  In  such  cases,  kinetic  parameters  were 
estimated  on  the  basis  of  R2  only.  Finally,  some  investigations 
also  report  kinetic  parameters  for  their  catalyst  formulation,  esti¬ 
mated  using  a  method  of  their  own  choosing.  For  each  such  case, 
the  parameter  values  obtained  using  our  technique  are  compared 
against  those  using  the  original  authors’  technique(s)  as  part  of  our 
discussion. 


Ci 


Qj+1  Jk 


QjJk 


55(e' 


QjJk 


+  IZi 


4j+\Jk 


)Sq 


Vie  {C0,C02,H2,H20,02} 
VjeZ,  0  <j  <  M-  1 
V/<  e  Z,  1  <  k  <  M 


(15) 


M  is  the  number  of  locations  within  the  reactor  for  which  con¬ 
centration  data  is  estimated.  Consequently,  8q  =  \/M=  qJ+i  -  qj, 
and  qo  =  0.  C\  is  the  experimentally  reported  inlet  concentration  of 
species  i  (i.e.  at  g  =  0). 

Our  overall  parameter  estimation  problem  is  concisely  stated  as 
follows: 

min  [Eq.  (12)]  subject  to  :  [Eq.  (15)],  (16) 

v 

This  is  a  highly  non-linear  (and  non-convex)  mathematical  pro¬ 
gram,  and  is  best  solved  using  numerical  optimization  techniques. 
For  this  investigation,  all  such  parameter  estimation  problems 
were  solved  using  the  IPOPT  optimization  algorithm  developed  by 
Wachter  and  Biegler  [15]. 

2.4.  Case-studies 

To  demonstrate  the  flexibility  and  utility  of  our  parameter 
estimation  technique,  we  applied  the  method  to  estimate  kinetic 


3.  Results  and  discussion 

In  Table  3,  we  report  values  of  R 2  and  R 2  obtained  by  fitting 
the  conversion  and  selectivity  measurements  reported  for  the  cat¬ 
alysts  in  each  of  the  eight  chosen  investigations  using  parameters 
estimated  through  our  approach  for  each  of  the  four  kinetic  mod¬ 
els  described  previously.  To  more  clearly  visualize  these  results, 
we  have  prepared  Fig.  2.  Here,  values  of  conversion  and  selectiv¬ 
ity  resulting  from  parameter  estimation  using  the  kinetic  model 
identified  in  Table  3  to  provide  the  highest  quality  fit  are  plotted 
against  the  conversion  and  selectivity  values  taken  from  the  origi¬ 
nal  studies.  The  identity  function  is  included  in  each  plot  as  an  ideal 
reference. 

In  many  cases  using  models  considering  WGS  activity,  it  was 
found  that  the  most  optimal  fit  of  the  observed  data  was  obtained 
when  the  WGS  reaction  was  in  fact  predicted  to  be  inactive 
(i.e.  ko  was  estimated  to  be  zero,  or  ED  was  estimated  to  be 
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Table  2 

Summary  of  the  eight  recent  literature  investigations  considered  in  this  study. 


Reference 

Catalyst 

Feed  composition3 

M 

Q 

Selectivity? 

Parameters? 

[4] 

1%  Pt  on  mordenite 

1%C0,  20%C02,  40%H2, 1.5%02 

11 

1 0,000  h-1 

No 

Yes 

[2] 

0.5%  Pt/0.02%  Fe,  on  alumina 

1%C0,  22%C02,  72%H2, 1%02 

8 

260  b-1 

Yes 

Yes 

[20] 

0.38%  Pt/10.56%  Mn,  on  alumina 

1%C0,  60%H2, 1%02 

6 

12,000  b”1 

Yes 

Yes 

[23] 

0.50%  Pt  on  alumina 

1%C0,  20%C02,  50%H2,  20%H20, 1%02 

11 

1 0,000  h”1 

Yes 

No 

[17] 

15%  CuO  on  ceria 

0.03%CO,  1%H2,0.03%02 

6 

83,000  b-1 

Yes 

Yes 

[18] 

5%  CuO  on  ceria/zirconia 

0.49%CO,  23.26%C02,  74.1 7%H2,  0.61  %02 

5 

SOOOh”1 

No 

No 

[21] 

0.5%  Ru  on  alumina 

0.5%CO,  18%C02,  37%H2,  5%H20,  0.5%02 

5 

67,000  b-1 

Yes 

No 

[19] 

2%  Au  on  thoria 

1%C0,  48%H2, 1%02 

7 

6000  b-1 

No 

No 

a  The  balance  of  the  indicated  feed  consists  of  inert  materials. 

Table  3 

Coefficients  of  determination  ( R 2 ,  R2)  listed  by  study  and  kinetic  model.  Best  fit  (largest  R2X  +  R2)  for  each  study  is  highlighted  in  boldface. 

Reference 

Kinetic  model 

Empirical 

Empirical  +  WGS 

Mechanistic 

Mechanistic  +  WGS 

[4] 

(0.986a) 

_b 

(0.9963) 

_b 

[2] 

(0.962,  0949) 

(0.970,  0.956) 

(0.904,  0.893) 

_b 

[20] 

(0.743,  0.753) 

_b 

(0.519,  0.860) 

_b 

[23] 

(0.921,  0.697) 

_b 

(0.663,  0876) 

_b 

[17] 

(0.981,0.911) 

_b 

(0.978,  0.903) 

_b 

[18] 

(0.887a) 

(0.9723) 

(0.9723) 

_b 

[21] 

(0781,  0.673) 

_b 

(0.597,  0.741) 

_b 

[19] 

(0.999a) 

_b 

(0998a) 

_b 

a  Because  no  selectivity  data  was  reported  in  the  referenced  study,  R 2  was  not  calculated  or  used  as  a  basis  for  parameter  estimation. 
b  Optimal  solution  has  zero  WGS  activity. 


prohibitively  large).  In  such  cases,  the  model  is  kinetically  iden¬ 
tical  to  the  associated  basic  model  (where  no  WGS  activity  was 
considered).  Consequently,  no  results  are  entered  for  such  cases  in 
Table  3. 

Table  4,  shows  the  kinetic  parameter  values  obtained  through 
estimation  for  each  (non-empty)  case  shown  in  Table  3.  Wherever 


possible,  parameter  values  reported  in  the  original  referenced  study 
are  included  beneath  our  estimated  values  for  comparison. 

From  our  results,  it  is  first  apparent  that  in  most  cases  the 
empirical  models  we  considered  provide  equivalent  or  better  fitting 
of  observed  catalyst  performance  than  do  the  mechanistic  mod¬ 
els.  In  fact,  in  only  one  case  ([4])  is  it  shown  that  a  mechanistic 


Table  4 

Estimated  kinetic  parameter  values  listed  by  study  and  kinetic  model.  Where  available,  comparison  parameters  for  a  given  study  are  included  in  italics.  Activation  energies 
(Ec,  Eh,  Ed,  El,\,  and  ELt 2)  are  reported  in  kj mol-1,  k  values  are  reported  in  the  following  units,  kc:  1/s  (mol  m-3)1  (a+b),  kH:  l/s  (mol  m-3)1  c,  kD:  1/s  (mol  m-3)  \  Icl,i : 
1  /s,  ki^ :  1  Is  (mol  m-3)  .  Exponents  (a,  b,  and  c)  are  unitless. 


Reference 

(4,r2) 

kc 

Ec 

l<H 

Eh 

a 

b 

c 

Kinetic  model:  empirical 

2*[4] 

210.986) 

5.17  xlO5 

31.41 

1.03  xlO7 

44.67 

0.78 

0.71 

1.20 

-a 

40.1 

-a 

-a 

-0.69 

0.68 

-a 

[2] 

(0962,0.949) 

38.95 

29.40 

1.23  xlO7 

55.87 

-0.78 

0.82 

3.00 

2*[20] 

210743,  0.753) 

8.76  x  106 

46.61 

2.85  xlO7 

55.44 

-1.31 

3.00 

0.95 

-a 

58.1 

-a 

101.8 

-a 

-a 

-a 

[23] 

(0.921,0.697) 

1.14  x  106 

48.44 

3.88  x  107 

56.97 

-0.73 

1.29 

1.20 

2*[17] 

210.981,0.911) 

1.88  xlO7 

29.86 

3.26  x  105 

60.41 

0 

2.02 

0 

-a 

48.5 

-a 

-a 

_a 

-a 

-a 

[18] 

(0.887) 

4.08  xlO5 

40.74 

1.34 

0.53 

0.97 

0 

0.52 

[21] 

(0.781,0.673) 

1.33  x  106 

35.56 

3.34  x  107 

49.10 

0.01 

1.07 

0.92 

[19] 

(0.999) 

4.90  x  104 

26.75 

4.93  x  105 

29.31 

-1.07 

1.34 

1.09 

Reference 

(R2X,R2) 

kc 

Ec 

kH 

Eh 

1<D 

Ed  o 

b 

c 

Kinetic  model:  empirical  +  WGS 

2*[2] 

210970,  0.756) 

0.246 

14.65 

758.8 

28.89 

1.86  xlO5 

81.93  -1.39 

0.79 

1.65 

352.8 

33.09 

20.53 

18.74 

4402 

34.10  -0.10 

0.50 

0.50 

[18] 

(0.972) 

1.93  xlO5  41.51 

3.99  x  10 

3  0.13 

1.58  xlO4 

40.79  0.29 

0.33 

2.97 

Reference 

(R^R2) 

kt,  i 

Ela 

kL,  2 

El,2 

kH 

Eh 

C 

Kinetic  model:  mechanistic 

[4] 

(0.996) 

5.72  x  104 

15.02 

4.16  x  10”3 

1.96  x  10”5 

1341.04 

14.59 

1.10 

[2] 

(0.904,  0.893) 

4.55  x  106 

59.64 

2.60  xlO7 

110.75 

2.23  xlO6 

54.48 

2.08 

[20] 

(0.519,  0.860) 

1.06  x  106 

37.43 

0 

0 

5.42  x  107 

57.68 

0.99 

[23] 

(0.663,  0.876) 

6.71  x  105 

37.91 

1.21  xlO4 

71.41 

2.68 

5.98 

0 

[17] 

(0.978,  0.903) 

1.95  x  108 

36.76 

0 

0 

9.21  x  104 

55.78 

0 

[18] 

(0.972) 

5.09  x  106 

29.62 

1.80  x  10”2 

1.90  x  10”7 

1.82  x  106 

49.48 

0.81 

[21] 

(0.597,  0.741) 

4.62  x  107 

39.73 

0 

0 

1.54  x  107 

46.72 

0.90 

[19] 

(0.998) 

3.36  x  106 

21.26 

8.19  xlO”3 

7.91  x  10”6 

5583.56 

25.45 

0.92 

a  No  value  for  the  indicated  parameter  was  provided  in  the  referenced  study. 
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Measured  Value 


Measured  Value 


Fig.  2.  Visualization  of  fitting  quality  for  each  chosen  study.  Conversion  (O)  and  selectivity  (x)  values  taken  from  the  original  study  are  plotted  along  the  horizontal  axis 
(“measured  value”);  values  obtained  through  parameter  estimation  are  plotted  along  the  vertical  axis  (“modeled  value”).  The  particular  kinetic  model  used,  and  the  resulting 
values  of  (R2,,  R2)  are  indicated  in  each  panel. 


model  provides  a  marginally  better  approximation  of  catalyst  per¬ 
formance. 

Mechanistic  models,  when  they  can  be  fit  to  observed  perfor¬ 
mance  data,  usually  generate  parameter  values  more  closely  tied 
to  the  underlying  chemistry,  and  thus  can  offer  more  detailed 
insight  than  can  general-purpose  empirical  models.  Often  however, 


such  detailed  information  is  not  needed  for  process  engineering 
applications  where  it  is  more  important  to  accurately  model 
macroscopic  catalyst  and  reactor  performance  than  it  is  to  accu¬ 
rately  model  microscopic  surface  interactions  between  reactants 
and  catalyst.  Moreover,  if  a  model  based  upon  an  inappropriate 
surface-reaction  mechanism  is  chosen,  it  may  be  impossible  to 
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accurately  model  the  catalyst  performance  and/or  the  parameter 
values  obtained  may  not  characterize  the  underlying  chemistry  as 
accurately  as  is  hoped. 

For  five  of  our  chosen  studies  ([4,2,17-19]),  the  mechanistic 
model  we  have  chosen  -  based  upon  a  bi-molecular  LH  surface 
mechanism  -  provides  a  good  approximation  of  the  overall  catalyst 
performance  (R%,  >  0.9,  R|  ^  0.9).  For  the  remaining  three  studies, 
however,  this  model  is  apparently  ill-suited  to  the  catalyst. 

We  also  note  that  in  two  cases  ( [20,2 1  ] ),  neither  model  (emprical 
or  mechanistic)  appears  capable  of  generating  parameters  that  very 
accurately  describe  the  observed  catalyst  performance  ( R2X  &  0.75, 
R 2S  &  0.7).  Although  the  empirical  model  still  offers  the  better  fit 
in  each  case,  the  results  suggest  that  a  better  solution  might  be 
achievable  using  a  kinetic  model  not  considered  here.  While  it  is 
generally  well  accepted  that  the  oxidation  of  CO  proceeds  via  an  LH- 
type  mechanism  over  Pt  [12],  the  catalyst  formulations  described 
in  these  two  studies  are  based  upon  Ru  [21]  and  a  Pt/MnOx  blend 
[20].  It  is  therefore  not  unreasonable  to  expect  that  their  surface 
chemistries  may  follow  different  models. 

In  general,  unless  a  precise  reaction  mechanism  is  known,  there 
appears  to  be  little  practical  advantage  in  modeling  the  catalytic 
oxidation  of  CO  using  mechanistic  models  (over  general  empirical 
models)  for  process  engineering  applications. 

The  results  in  Table  3  also  provide  some  insight  into  the  role  of 
the  WGS  reaction  within  PROX  reactors.  In  most  cases,  the  best  fit  of 
the  experimental  conversion  and  selectivity  data  is  obtained  when 
this  reaction  is  inactive.  In  these  cases,  therefore,  the  observed  con¬ 
version  and  selectivity  is  adquately  explained  through  catalysis  of 
the  principal  CO  and  H2  oxidation  reactions  only.  For  two  of  our  cho¬ 
sen  studies,  however,  the  best  fit  to  observed  catalyst  performance 
is  obtained  with  non-zero  WGS  activity. 

Choi  and  Stenger  [2]  reported  that  WGS  activity  was  expected 
to  be  most  significant  under  conditions  where  large  amounts  of  H2 
and  C02  are  present  in  the  feed,  overwhelming  the  forward  reaction 
as  favored  by  chemical  equilibrium.  It  is  therefore  not  surprising 
that  the  two  cases  for  which  non-zero  WGS  activity  was  predicted 
involve  feed  streams  containing  at  least  90%  H2/C02.  Our  results 
are  therefore  consistent  with  the  expectations  of  Choi  and  Stenger, 
and  demonstrate  that  simultaneous  catalysis  of  the  WGS  reaction 
must  be  considered  when  feeding  PROX  with  a  large  amounts  of  H2 
or  C02. 

Examining  the  parameter  values  given  in  Table  4,  we  first  find 
that,  in  general,  the  activation  barrier  for  CO  oxidation  (£c)  is  lower 
than  that  for  FI2  oxidation  (Eh).  This  is  to  be  expected  from  catalysts 
that  preferentially  oxidize  CO. 

Examining  the  parameter  data  further,  we  discover  that  trends 
in  parameter  values,  even  for  catalyst  preparations  based  upon  the 
same  active  metal  or  metal  oxide,  are  not  immediately  evident. 
Because  the  apparent  chemical  kinetics  of  a  catalyst  are  affected  not 
only  by  its  active  component,  but  also  its  support,  its  preparation 
method,  and  the  operating  conditions  to  which  it  is  exposed  [  1 6,22  ] 
this  lack  of  discernable  trends  is  not  discomforting. 

Finally,  we  note  that  three  of  our  chosen  studies  ([4,18,19])  did 
not  report  selectivity  data  as  a  function  of  temperature.  Kinetic 
parameters  in  these  cases  could  therefore  be  estimated  only  on 
the  basis  of  CO  conversion.  As  a  result,  the  importance  of  the  H2 


oxidation  kinetics  are  de-emphasized  in  these  cases,  and  we  there¬ 
fore  have  less  confidence  in  the  accuracy  of  the  kinetic  parameters 
I<h,  Eh,  and  C.  This  is  particularly  true  in  the  case  of  Ratnasamy 
et  al.  [18]  for  which  the  emipirical  model  estimates  a  value  of  EH 
substantially  less  than  Ec;  values  inconsistent  with  the  preferential 
oxidation  of  CO. 
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